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ABSTRACT 

This  work  derives  the  Cramer-Rao  lower  bound  (CRLB)  for  an  acoustic  target  and  sensor  localization  system 
in  which  the  noise  characteristics  depend  on  the  location  of  the  source.  The  system  itself  has  been  previously 
examined,  but  without  deriving  the  CRLB  and  showing  the  statistical  efficiency  of  the  estimator  used.  Two 
different  versions  of  the  CRLB  are  derived,  one  in  which  direction  of  arrival  (DOA)  and  range  measurements 
are  available  (“full-position  CRLB”),  and  one  in  which  only  DOA  measurements  are  available  (“bearing-only 
CRLB”).  In  both  cases,  the  estimator  is  found  to  be  statistically  efficient;  but,  depending  on  the  sensor-target 
geometry,  the  range  measurements  may  or  may  not  significantly  contribute  to  the  accuracy  of  target  localization. 

Keywords:  CRLB,  FIM,  estimator  efficiency,  acoustic  localization,  fusion 

1.  INTRODUCTION 

In  any  estimation  system  the  ultimate  goal  is  to  extract  the  maximum  information  from  the  available  data.  The 
Fisher  information  matrix  (FIM)  provides  a  measure  of  the  total  information  available  from  the  observations  of 
the  system,  and  its  inverse  provides  the  Cramer-Rao  lower  bound  (CRLB)  pQ.  A  statistically  efficient  estimator 
is  one  in  which  the  variance  of  the  estimation  error  meets  the  CRLB,  and,  therefore,  extracts  all  of  the  available 
information  from  the  observations. 

The  CRLB  and  statistical  efficiency  of  an  acoustic  localization  system  will  be  examined  here,  based  on  the 
system  described  in  [BJ,  which  is  meant  to  estimate  the  location  of  the  source  of  a  detected  gunshot.  Each  sensor 
node  of  the  system  is  assumed  to  provide  an  estimated  bearing  (DOA)  to  the  target,  and,  if  the  sensor  node  lies 
within  a  certain  field-of-view  (FOV),  a  range  estimate  and  bullet  trajectory  estimate  as  well.  For  those  sensors 
which  provide  estimated  range,  the  noise  term  will  be  dependent  on  the  position  of  the  source.  Each  sensor 
node’s  estimates  are  passed  to  a  fusion  center  to  perform  the  overall  estimation  of  the  target  position  by  fusing 
all  the  local  estimates. 

A  number  of  papers  have  examined  the  problem  of  target  localization  in  passive  sensor  environments,  including 
[2M1I51I51I10H12|.  The  work  of  m  generalizes  the  results  of  Rj  to  include  sensor  position  uncertainty,  however, 
neither  paper  examines  the  CRLB  to  see  whether  the  estimator  is  statistically  efficient.  In  [21151151112] .  different 
applications  of  localization  with  passive  sensors  are  studied  which  also  consider  the  CRLB.  However,  in  [T2], 
no  estimation  scheme  is  shown  to  meet  the  CRLB.  In  [5J  the  ML  estimation  scheme  examined  is  shown  to  be 
statistically  efficient  only  when  a  significant  number  of  measurements  are  utilized.  In  none  of  the  previously 
mentioned  papers  were  cases  of  position  dependent  measurement  noise  considered. 

In  this  work,  the  CRLB  of  the  central  estimator  (fuser)  is  derived  for  two  cases:  a  “bearing-only”  case,  which 
assumes  that  no  range  estimates  are  available  from  the  sensor  nodes,  and  a  “full-position”  case,  which  assumes 
that  range  estimates  are  available  in  addition  to  the  DOA. 

Section  [2]  provides  an  overview  of  the  system  in  question.  Section  [3]  provides  the  expressions  necessary  to 
evaluate  the  CRLB  for  the  problem  in  question,  both  with  and  without  the  position  dependent  noise  terms. 
Section  [4]  describes  the  simulation  scenarios  and  provides  the  results.  Finally,  Section  [5]  concludes  the  paper. 
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2.  LOCALIZATION  SYSTEM  OVERVIEW 


The  system  to  be  examined  here  will  be  identical  to  the  one  described  in  {6] .  A  brief  overview  of  the  system  will 
be  provided  here  to  introduce  the  concepts  and  notations. 

A  number  of  acoustic  sensors  are  placed  throughout  a  surveillance  region  with  the  intent  of  detecting  gunfire 
and  estimating  the  position  of  the  source.  The  target  (source)  location  will  be  denoted  as 


T  = 


Tx 
Ty . 


and  the  ith  sensor  location  will  be  denoted  as 


The  problem  will  be  assumed  constrained  to  a  two-dimensional  plane  for  simplicity. 


(1) 

(2) 


2.1  Sensor  Nodes 

Each  sensor  will  provide  at  most  five  measurements  (local  estimates):  the  DO  A  angle  to  the  shooter  pi,  based 
on  the  detection  of  the  muzzle  blast,  the  DOA  angle  of  the  shockwave  from  the  bullet  pi,  the  time  difference  of 
arrival  (TDOA)  between  the  muzzle  blast  and  the  shockwave  A,  and  the  sensor  location  Six  and  Si  (through 
a  GPS  sensor  at  each  node).  Each  measurement  is  assumed  to  be  corrupted  by  zero-mean  Gaussian  noise, 
with  standard  deviations  of  a^,  a^,  aT ,  OiT  and  aly  ,  respectively.  The  shockwave  (and  therefore  the  TDOA 
measurement)  will  only  be  visible  to  sensor  nodes  which  are  within  a  limited  field-of-view  (FOV)  around  the 
path  of  the  bullet.  The  FOV  will  be  tt  —  2d ,  where 

0  =  sin-1  (3) 


and  m  is  the  Mach  number  of  the  bullet,  assumed  here  to  be  m  =  2  [7|. 
The  target  bearing  from  the  ?’th  sensor  node  is 

'Ty~Siy 


Pi  =  tan 


-l 


Tx  -  S, 


and  the  DOA  angle  of  the  shockwave  is 


= 


—  f  —  9+UJ  if7T  +  CLl< 


T  Q  P)  if  -p  -p  Cc?  <C  7T  — (—  Cc? 


where  oj  is  the  angle  of  the  trajectory  (counter-clockwise)  with  respect  to  the  a;-axis. 
Sensor  i  will  pass  on  the  local  estimate^ 


Zi  = 

Pi 

and  Si  = 

<3- 

tUi_ 

(4) 

(5) 


(6) 


to  a  fusion  center,  where  pi  is  the  estimated  target  bearing,  f,  is  the  estimated  target  range  obtained  from  pi 
and  A  using  0,  and  Cbi  is  the  estimated  bullet  trajectory  angle.  In  view  of  0,  the  estimated  bullet  trajectory 
can  be  obtained  directly  from  pi  and  the  standard  deviation  of  o)j  will  be  av. 

The  estimated  range  will  be  [6] 

n  =  i - 71 — M  (7) 

1  -  cos(0i  -  (Pi) 


JWe  use  similar  notation  and  terminology  as  in  [6]. 
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where  c  is  the  speed  of  sound  (assumed  to  be  known  perfectly).  The  variance  of  the  range  estimate  —  which  is 
dependent  on  both  range  and  bearing  —  can  be  approximated  as  [6] 


o$(T,Si,w) 


c2&t  +  rf(cr2  +  sin2 (</>$  -  y Oj) 
(l  -  cos {(t>i  - 


(8) 


The  likelihood  functions  corresponding  to  the  estimates  zt  and  Si  are 

p(zi\T,Si,uj)  =  J\f(zi;/j,Zi,T,Zi) 


(9) 


where 


and 


where 


0 

0  ' 

= 

Ti 

and  YiZi  = 

0 

a*.(T,Si,u)) 

0 

UJ 

_  0 

0 

< 

p(Si\Si)=Af(Si-,Si,'ESi) 


(10) 

(11) 

(12) 


The  errors  in  ([9|)  and  m  are  assumed  to  be  uncorrelated  at  each  sensor  and  across  the  sensors. 

The  sensors  are  assumed  to  obtain  their  locations,  albeit  imperfectly,  from  GPS.  d  If  desired,  the  algorithm 
can  perform  simultaneous  localization  of  the  target  and  sensors.  For  the  sensors,  the  GPS  localization  serves  as 
a  prior  and  guarantees  complete  observability  for  the  target-sensor  complex.  The  final  estimates  of  the  sensor 
locations  can  be  slightly  improved  over  their  initial  GPS  estimates,  but  the  improvement  this  makes  to  the  final 
target  localization  is  negligible. 


2.2  Centralized  Fusion 

The  estimates  z%  and  S)  from  each  sensor  are  passed  on  to  a  fusion  center  in  order  to  determine  the  estimate  x 
by  means  of  the  Iterated  Least  Squares  (ILS)  estimator  [T.  The  parameter  vector  to  be  estimated  is 

n=[Tx  Ty  UJ  Slx  Sly  ...  Snx  S„y]’  (13) 


with  observations 


where 


(14) 


(15) 


In  order  to  utilize  the  ILS  estimation  algorithm,  an  initial  estimate  of  x  is  needed.  It  has  been  noted  [6]  that 
the  ILS  estimator  is  sensitive  to  the  initial  estimate  and  may  diverge  if  the  initial  estimate  is  too  far  from  the 
truth. 

While  the  initialization  of  the  target  position  could  be  performed  by  using  the  bearing  and  range  measure¬ 
ments  (using  the  nodes  with  range  measurements),  the  large  variance  of  the  range  measurements  was  found  to 
occasionally  cause  divergence  in  the  ILS  algorithm.  A  more  robust  initialization  was  found  to  follow  a  similar 

2If  the  sensor  position  estimates  contain  a  common  bias  across  sensors,  the  relative  sensor  registration  will  be  unaffected 
and  the  target  estimate  will  exhibit  the  same  bias. 
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method  to  that  used  in  [3].  This  method  of  initialization  will  utilize  only  the  available  bearing  measurements 
from  each  sensor  ©,  which  can  be  rewritten  as 


tan  —  1 

Slx  tan  0i  -  Si/ 

tan  02  —  1 

T  = 

S2x  tan  02  -  S2y 

tan  0n  —  1 

S nx  tan  0n  S ny 

and  T  can  be  solved  as 


T  =  Afb 


where  is  the  (right)  pseudo-inverse  of  A. 

Also,  note  that  (USD  can  be  rewritten  utilizing  the  expression 


<j)i  =  cot  1 


(16) 


(17) 


(18) 


which  is  simply  (j4])  rewritten  using  the  cotangent  function.  As  suggested  in  [3],  use  of  the  cotangent  function 
has  been  made  when  the  measured  bearing  is  between  45°  and  135°  or  between  -45°  and  -135°. 

To  complete  the  initialization  of  x,  Si  can  be  taken  as  Si  and  u  can  be  taken  as  the  average  of  cD.;. 


3.  CRAMER- RAO  LOWER  BOUND 


The  CRLB  provides  a  lower  bound  on  the  covariance  matrix  of  the  estimate  x  as 

E  [(x  —  x)  (x  —  x)']  >  J~l  (19) 

where  J  is  the  Fisher  information  matrix  (FIM) 

J  =  £{[VxA(x)][VxA(x)]'}  (20) 

and  A(x)  is  the  negative  log-likelihood  function  (NLLF) 

A(x)  =  ^2  2  (Zi  ~  fe(x)),S^(x)~1  (zi  ~  Mzi(x))  +  y  (Si  -  S»)  ^(^-S^+ln^/x))  (21) 

»= l  L 

with  the  unnecessary  constant  terms  omitted. 

In  view  of  ©,  (1101).  and  CEB).  the  NLLF  can  be  written  a3 


aw  =  oE 


2  l  r*  r'l 


l 


a  %  (T,  Si,  u>) 


(fi  -  nf  +  (w»  -  w)2  +  ( Six  -  Si})  +  [Siy  -  Si y) 


+  ln  {crl.(T,Si,w)) 


(22) 


Note  that  the  noise  term  for  the  ith  sensor  node’s  range  estimate  is  dependent  on  the  target-sensor  geometry. 

In  order  to  deal  with  the  position  dependent  noise  term,  the  “full-position”  FIM  corresponding  to  E^l)  will 
be  broken  into  “bearing-onlv  □  and  “range-only”  FIMs 

J  =  Jb  T  Jr  (2f>) 

where  Jb  denotes  the  “bearing-only”  FIM  calculated  when  ignoring  the  range  term  of  ((22]) :  and  Jr  denotes  the 
“range-only”  FIM  calculated  for  only  the  range  portion  of  (l22l) . 

Two  versions  of  the  CRLB  will  be  evaluated  in  the  sequel:  the  bearing-only  CRLB,  ,/ft~  ,  and  the  full-position 
CRLB,  J”1. 

3Note  that  uj  is  decoupled  from  the  remaining  terms  of  x,  and  the  ML  estimation  of  uj  amounts  to  simple  averaging, 
with  a  CRLB  of  crj/n. 

4This  is  a  slight  misnomer  since  the  sensor  position  terms  will  also  be  included  in  the  bearing-only  FIM. 
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3.1  Bearing-Only  FIM 

The  bearing-only  FIM  J ],  will  be  straightforward,  as  none  of  the  variances  in  this  case  are  dependent  on  x. 
The  gradient  of  Kill),  necessary  for  calculation  of  the  FIM,  is  detailed  in  Appendix  [B]  The  necessary  terms  to 
calculate  the  FIM  are 


f  d X  dX\ 
\dTxdTx) 


f  dX 

dX  ) 

{  dTx 

dTyj 

f  dX 

dX] 

l  &TX 

du  J 

f  dX 

dX  I 

l  dTx 

dSixj 

f  dX 

dX  \ 

dTx  8Sly  J 


f  d\  d\  \ 
{dfyWyj 


f  d\  d\  \ 

f  dX  OX  \ 
{ dTvdSuJ 


dX  dX  \ 
dTydS~J 
dXdX\ 

du>  8uj  J 
dX  dX  \ 
doj  dSix  J 
dX  dX  \ 
du;  dSiy  j 


E 

E 


dX  dX  \ 
dSlx  dSlx  J 
dX  dX  \ 
dS~dS~] 


f  dX  dX  \ 

{  ds~as~! 


The  matrix  Jt  can  then  be  constructed  with  the  appropriate  entries  from  (1241)  (l38l). 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


3.2  Range-Only  FIM 

Since  the  variance  of  the  range  estimate  from  each  sensor  node  (for  which  the  shooter  is  in  the  FOV  necessary 
to  obtain  a  range  estimate)  is  dependent  on  x,  the  expressions  for  the  entries  of  the  FIM  in  the  range-only  case 
are  more  complicated.  The  terms  for  the  gradient  of  (1221)  (including  the  range  terms)  are  detailed  in  Appendix 
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[Cl  The  necessary  terms  to  calculate  the  FIM  Jr  are 


E 

E 

E 

E- 

E- 

E 

E 

E- 

E- 


dX  d\ 


l  dTx 

dTx 

f  dX 

dX 

l  OT, 

dTy 

r  ox 

dX' 

{ m 

;  OOJ 

f  dX 

dX 

{0TX 

dSlx 

1  dX 

dX 

l  dTx 

dSly 

r  ox 

dX 

\dTy 

dTy 

r  ox 

dX' 

{or. 

,  du  _ 

r  dx 

dX 

l dTy 

dSlx 

r  ox 

dX 

l dTy 

dSiy 

dXdX\ 

du  du  J 


E 

E 

E< 

E< 

E< 


'  dX  dX  \ 

{ du  dSix  J 
'  dX  dX  \ 
l  du  dSiy  j 
dX  dX  \ 
dS~dS~J 
dX  dX  \ 

dsZdS~) 

dX  dX  \ 
dSly  dSly  ) 


- 2  X+2 + 51 


(Tx-Sixf 


r2cr2 

i=  1  i=  1  i  ri 


=  2  aibi  +  ^2 

i= 1  i= 1 

n 

=  2  ajCj 


(Tx  -  5,J(Ty  -  Sly) 


r2cr2 

l  Vi 


2q2  {Tx  —  Six)2 


r2<j2 

'  l  Vi 


—  2dibi 


(Tx  -  S*.)^  -  &  ) 


r^/T"2 


=  ?£«+E 

i=l  i= 1 

n 

=  2  y  '  6jCi 


i=l 

2(ijbi 

—26?  - 
n 

*E<* 

i= 1 
2  CliCi 


(Tx  -  Six)(Ty  -  Sly) 


r2a2 

l  Vi 


2 

- 9“^ - 

nr**  rr* 


=  -2b; 


Ci 


=  2  af 


2  ,  (T*-^)2 


r2(j2 

'  i  ri 


2 &ibi 
26?  + 


(T,  -  5,J(Ty  -  &  ) 


r2a2 

'  l  V  i 


C Ty-Siy )2 


7*?(j2 
'  %  w  r. 


E 

E 

E 


dX  dX 
dS~dS~ 

dX  dX 
OS~dS~ 

dX  dX 

ds~ds~ 


=  0 

=  0 

ilXj 

=  0 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 
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where 


a  _  sin (<fc  -  ipi)  /sign (Tx  -  Six){al  +  a^n  sin(ygi)  _  Ty  -  Siy  \ 

1  -  cos  {(pi  -  ipi)  V  (1  -  cos  {fa  -  .  rf  y 

_  sin (fc  -  <pi)  /sign (Tx  -  5)J(a|  +  cos(^)  _  -  Sia!  \ 

i  -  cos (</>*  -  v?;)  V  (i  -  cos(y-i  -  <^))o-2.  rf  y 

_  rfjal  +  (7y )  sin(<fo  -  ipi)  cosjfa  -  <^)  _  sin(^  -  w) 

(1  -  cos (4>i  -  ipi))2^.  1  -  cos (</>j  -  (^i) 

The  matrix  Jr  can  then  be  constructed  with  the  appropriate  entries  from  (13911  (1561). 


(57) 

(58) 

(59) 


4.  SIMULATION  RESULTS 

The  simulation  scenarios  examined  here  include  the  scenarios  of  [6]  and  an  additional  modified  scenario  with 
fewer  sensors.  For  each  scenario,  the  Mach  number  of  the  bullet  is  assumed  to  be  m  =  2,  and  the  speed  of  sound 

is  assumed  to  be  c  =  342m/sec.  The  measurement  noise  standard  deviations  are  =  4°,  ay  =  1ms,  and 

(7 ix  =  <Jiy  =  2  m.  The  simulations  were  performed  for  100  Monte  Carlo  runs  for  each  scenario. 

In  Scenarios  1  and  2,  there  are  five  sensor  nodes  located  at 

[127  20  90  136  182]  ,  . 

*  -  107  22  0  68  59  ('bU^ 


In  Scenario  1,  the  target  is  located  at  T  =  [50  50]/  and  the  bullet  is  fired  at  a  trajectory  of  uj  =  30° 

(counter-clockwise  from  the  cc-axis).  Due  to  the  location  of  the  sensors  and  the  trajectory  of  the  bullet,  only 
sensors  1,  4  and  5  receive  the  shockwave  and  are  able  to  send  range  and  bullet  trajectory  estimates  to  the  fusion 
center. 

The  results  of  Scenario  1  are  shown  in  Figure  [T|  and  [2]  Each  figure  shows  the  true  locations  of  the  target 
and  sensors,  along  with  the  corresponding  95%  error  ellipses.  Figure  Q]  shows  the  error  ellipses  corresponding  to 
the  sample  covariance  matrix  calculated  from  the  estimation  errors  over  the  100  Monte  Carlo  runs  when  range 
measurements  are  not  available  at  the  fusion  center  (dashed  line),  and  the  covariance  matrix  from  the  bearing- 
only  CRLB  (solid  line).  Figure [2] shows  the  covariance  matrix  calculated  from  the  estimation  errors  when  both 
range  and  bearing  data  are  available  at  the  fusion  center,  and  the  covariance  matrix  from  the  full-position  CRLB. 

The  covariance  matrices  from  the  bearing-only  CRLB  and  the  full-position  CRLB  closely  match  the  covari¬ 
ances  of  the  estimation  errors  calculated  from  the  simulation.  This  first  indicates  that  the  ILS  estimation  carried 
out  by  the  fusion  center  is  statistically  efficient.  Additionally,  the  fact  that  the  two  CRLB  matrices  closely  match 
suggests  that  very  little  information  is  gained  from  the  range  estimates  sent  from  sensors  1,  4  and  5.  This  is  due 
to  the  fact  that  the  target-sensor  geometry  is  such  that  the  target  bearing  measurements  alone  can  achieve  good 
localization. 

In  Scenario  2,  the  target  is  located  at  T  =  [150  — 50]  ^  and  the  bullet  is  fired  at  a  trajectory  of  w  =  170°. 

Due  to  the  location  of  the  sensors  and  the  trajectory  of  the  bullet,  only  sensors  2  and  3  receive  the  shockwave 
and  are  able  to  send  range  and  bullet  trajectory  estimates  to  the  fusion  center. 

The  results  of  Scenario  2  are  shown  in  Figures  [3]  and  [4]  Each  figure  once  again  shows  the  various  95% 
error  ellipses.  Figure  [3]  shows  the  error  ellipses  of  the  bearing-only  CRLB  and  the  estimation  errors  when  range 
measurements  are  not  available  at  the  fusion  center.  Figure  2]  shows  the  error  ellipses  of  the  full-position  CRLB 
and  the  estimation  errors  when  both  range  and  bearing  data  are  available  at  the  fusion  center. 

The  covariance  matrices  from  both  cases  of  the  CRLB  again  closely  match  those  obtained  from  the  estimation 
errors,  indicating  that  the  estimator  is  once  again  efficient  and  the  range  estimates  carry  very  little  information. 

Scenario  3  consists  of  an  identical  situation  to  Scenario  1,  but  with  sensors  2  and  3  removed.  In  this  case, 
the  geometry  of  the  sensors  and  target  are  poor,  with  each  sensor  having  very  similar  line-of-sight  (LOS)  angles 
to  the  target. 
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Figure  1:  Scenario  1,  bearing-only  CRLB  ellipses  and 

error  ellipses  of  estimated  target  and  sensor  locations 
(all  95%). 
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Figure  2:  Scenario  1,  full-position  CRLB  ellipses  and 

error  ellipses  of  estimated  target  and  sensor  locations 
(all  95%) 


Figure  3:  Scenario  2,  bearing-only  CRLB  ellipses  and 
error  ellipses  of  estimated  target  and  sensor  locations 
(all  95%). 


Figure  4:  Scenario  2,  full-position  CRLB  ellipses  and 
error  ellipses  of  estimated  target  and  sensor  locations 
(all  95%). 


Figure  5:  Scenario  3,  bearing-only  CRLB  ellipses  and 
error  ellipses  of  estimated  target  and  sensor  locations 
(all  95%). 
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Figure  6:  Scenario  3,  full-position  CRLB  ellipses  and 
error  ellipses  of  estimated  target  and  sensor  locations 
(all  95%). 


Table  1:  Normalized  Estimation  Error  Squared  (NEES)  (95%  probability  region  is  [1.63,  2.41]). 


NEES  (Tbo) 

NEES(f) 

NEES(f)  wrt  CRLBbo 

Scenario  1 

1.95 

2.10 

1.84 

Scenario  2 

2.16 

1.97 

2.00 

Scenario  3 

2.12 

2.20 

1.41 

The  results  of  Scenario  3  are  shown  in  Figures  [5]  and  [6]  In  this  case,  the  error  ellipses  for  the  bearing-only 
CRLB  and  the  full-position  CRLB  no  longer  coincide.  The  covariance  of  the  estimation  errors  closely  match 
their  respective  CRLB  matrices,  again  indicating  that  the  estimator  is  efficient. 

Additionally,  the  normalized  estimation  error  squared  (NEES)  for  the  source  localization  was  examined  for 
each  scenario,  using  both  the  bearing-only  CRLB  and  the  full-position  CRLB  to  provide  a  statistical  confirmation 
of  the  efficiency  of  the  estimator.  In  each  case,  the  CRLB  was  evaluated  at  the  true  x.  The  NEES  was  calculated 
for 

(i)  the  fused  position  estimation  errors  using  bearing-only  measurements  with  the  bearing-only  CRLB, 

(ii)  the  fused  position  estimation  errors  using  full-position  measurements  with  the  full-position  CRLB,  and 

(iii)  the  fused  position  estimation  errors  using  full-position  measurements  with  the  bearing-only  CRLB  (to 
quantify  the  benefit  of  range  measurements). 


The  NEES  results  (with  the  95%  probability  region  based  on  the  chi-square  distribution  with  two  degrees  of 
freedom  and  100  Monte  Carlo  runs  being  [1.63,  2.41])  are  in  Table  [T]  The  first  two  columns  (case  [(I)]  and 


(ii)  above)  show  that  in  all  scenarios,  with  bearing-only  or  full  position,  the  fused  estimate  of  the  target  (source) 


position  matches  the  CRLB,  i.e.,  it  is  statistically  efficient. 


The  third  column  (case  (iii)  above)  shows  that  the  target  position  errors  with  full  position  measurements 
match  in  Scenario  1  and  2  the  bearing-only  CRLB,  i.e.,  the  range  availability  does  not  contribute  to  the  accuracy 
of  the  fused  estimate.  However,  in  Scenario  3,  the  availability  of  the  range  reduces  the  MSE  by  30%  compared 
to  the  bearing-only  case. 


Additionally,  the  aforementioned  scenarios  were  simulated  for  multiple  levels  of  angular  measurements  noise. 
The  standard  deviations  a ^  and  were  varied  from  10%  to  150%  of  their  original  value  of  4°.  The  remaining 
noise  standard  deviations  remained  the  same  as  in  the  previous  simulations. 
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Figure  7:  Scenario  1,  NEES  for  different  levels  of  angular  Figure  8:  Scenario  2,  NEES  for  different  levels  of  angular 
measurement  noise.  measurement  noise. 


Measurement  Noise  Standard  Deviation  Scaling  Factor 


Figure  9:  Scenario  3,  NEES  for  different  levels  of  angular  measurement  noise. 


Figure  [7]  shows  the  NEES  of  Scenario  1  (from  Table  [U  along  with  the  95%  confidence  regions.  Figures  [8]  and 
[9]  show  identical  figures  for  Scenarios  2  and  3,  respectively. 

Figures  [7]  [9]  mirror  the  results  of  Tabled  demonstrating  statistical  efficiency  of  the  fused  target  position 
estimate,  with  the  exception  of  very  inaccurate  angular  measurements  and  poor  sensor-target  geometry. 

Figure  [10]  shows  the  target  position  RMSE  (coordinate-combined)  for  Scenario  1,  along  with  the  CRLB,  for 
both  the  bearing-only  and  full-position  cases.  This  figure  also  shows  the  95%  confidence  region  around  the  CRLB 
for  the  sample  RMSE  over  100  Monte  Carlo  runs.  Figures  ITT1  and  ITT!  show  identical  figures  for  Scenarios  2  and 
3,  respectively. 

Figures  [TUI  and  fill  show  that,  over  a  range  of  angular  measurement  noise  levels,  the  favorable  geometry  of 
Scenarios  1  and  2  provide  for  very  little  differentiation  in  the  performance  of  target  localization  with  or  without 
range  measurements.  Figure IT2l  shows  that,  for  the  less  favorable  geometry  of  Scenario  3,  the  inclusion  of  range 
measurements  provides  a  significant  increase  in  the  accuracy  of  target  localization. 

Additionally,  the  two  versions  of  the  CRLB  (bearing-only  and  full-position)  can  be  compared  to  gain  insights 
into  a  particular  scenario.  The  ratio  of  the  area  of  the  bearing-only  CRLB  ellipse  to  the  full-position  CRLB  ellipse 
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RMSE 


Figure  10:  Scenario  1,  target  position  RMSE  for  different  levels  of  angular  measurement  noise  and  95%  probability 
region  around  the  CRLB. 


Measurement  Noise  Standard  Deviation  Scaling  Factor 
Figure  11:  Scenario  2,  target  position  RMSE  for  differ¬ 
ent  levels  of  angular  measurement  noise  and  95%  prob¬ 
ability  region  around  the  CRLB. 


Measurement  Noise  Standard  Deviation  Scaling  Factor 
Figure  12:  Scenario  3,  target  position  RMSE  for  differ¬ 
ent  levels  of  angular  measurement  noise  and  95%  prob¬ 
ability  region  around  the  CRLB. 
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Figure  13:  Comparison  of  CRLB  and  CRLBbo  (|  J^1  J|1,/2)  over  a  2-d  grid  of  shooter  position  (for  a  fixed  aimpoint 
at  *). 
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can  be  calculated  as  |  This  is  plotted  over  a  two-dimensional  grid  corresponding  to  various  shooter 

locations  in  Figure[13l  for  sensor  locations  identical  to  Scenario  3.  Figure [T3] clearly  shows  the  shooter  locations 
where  the  range  measurements  are  most  beneficial  and  the  bearing-only  localization  will  perform  particularly 
poorly. 


5.  CONCLUSIONS 

The  CRLB  and  statistical  efficiency  were  examined  for  multiple  scenarios  of  a  localization  system  with  position 
dependent  noise  terms;  specifically,  a  recently  developed  system  utilizing  acoustic  gunfire  detection  sensors  jjl]. 
The  CRLB  was  derived  for  both  the  case  of  sensor  nodes  which  send  both  bearing  and  range  measurements  ( “full- 
position”  CRLB)  and  those  that  do  not  (“bearing-only”  CRLB).  In  cases  where  the  sensor-target  geometry 
is  favorable  for  angle-only  localization,  the  bearing-only  CRLB  closely  matches  the  actual  estimation  errors, 
suggesting  that  there  is  little,  if  any,  information  contained  in  the  range  measurements  in  those  cases.  If  the 
geometry  is  poor,  however,  as  in  Scenario  3,  the  full-position  CRLB  is  the  only  one  that  will  closely  match  the 
actual  estimation  errors.  The  results  show  both  that  the  estimator  used  in  this  particular  acoustic  localization 
system  is  efficient,  and  that  the  full-position  CRLB  could  act  as  an  accurate  means  of  performance  prediction 
for  such  a  system. 

APPENDIX  A.  GENERAL  FIM  FOR  MULTIVARIATE  GAUSSIAN  CASE 

The  FIM  for  the  case  of  a  multivariate  Gaussian  likelihood  is  as  follows  [9] . 

Assuming  the  likelihood  of  the  parameter  vector  x  given  observations  y  is  distributed  as  A /”(/r(x),  £(x)),  the 
NLLF  (ignoring  the  irrelevant  constant  terms)  is 

A(x)  =  iln|E(x)|  +  i(y-^(x))'S(x)-1(y-^(x))  (61) 
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The  gradient  terms  of  the  FIM  are 


2M  =  hrfE-VaE(x)i  1 
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_  dxi 

where 

The  (i,  j)th  entry  in  the  FIM  is  then 


X  =  [xi  X2 
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1  ’  dxn 
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(64) 


APPENDIX  B.  GRADIENT  TERMS  FOR  BEARING-ONLY  FIM 

The  terms  of  the  gradient  of  the  log-likelihood  function  (1221).  when  ignoring  the  range  measurements,  are 
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APPENDIX  C.  GRADIENT  TERMS  FOR  FULL-POSITION  FIM 

The  terms  of  the  gradient  of  the  log-likelihood  function  (1221).  including  the  range  measurements,  are 
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where 
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